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ABSTRACT 

The topology of weak lensing fields is studied using the 2-dimensional genus statistic. Simulated fields 
of the weak lensing convergence are used to focus on the effect of nonlinear gravitional evolution and to 
model the statistical errors expected in observational surveys. For large smoothing angles, the topology 
is in agreement with the predictions from linear theory. On smoothing angles smaller than 10', the genus 
curve shows the non-Gaussian signatures of gravitational clustering and differs for open and flat cold 
dark matter models. Forthcoming surveys with areas larger than 10 square degrees should have adequate 
signal-to-noise to measure the non-Gaussian shape and the Sl-dependence of the genus statistic. 

Subject headings: cosmology: theory — gravitational lensing — large-scale structure of universe — 
methods: numerical 



1. INTRODUCTION 

Mapping the mass density field with weak gravitational 
lensing is a promising method to analyze large-scale struc- 
ture in the universe. The statistical nature of the mass 
density field is directly predicted by theories of the early 
universe. For instance, most infiationary theories predict 
a random Gaussian initial mass density field, while struc- 
ture formation by cosmic defects predicts a highly non- 
Gaussian field. As the density field evolves, gravitational 
clustering produces characteristic non-Gaussian features. 
Thus measuring the non-Gaussianity of the density field 
can be an important test of the origin and evolution of 
large-scale structure. 

The genus statistic is known to be a useful topological 
discriptor of the the clustering pattern of the mass density 
field (Gott, Melott & Dickinson 1986; Gott et al. 1989). 
The "genus curve", which plots the measured genus ver- 
sus the threshold for iso-density contours, has a univer- 
sal shape for a random Gaussian field. Departures from 
this Gaussian shape reveal the mechanism of the forma- 
tion of nonlinear structures. So far, the genus statistic has 
been applied to the distribution of galaxies, and of clus- 
ters (Moore et al. 1992; Rhoads, Gott & Postman 1994; 
Canavezes et al. 1998). Since these tracers are affected by 
biasing, it is desirable to directly use the mass density field 
for the genus statistic. Fortunately, recent developments 
in measuring weak lensing in wide area surveys suggest 
that this may be feasible in the near future (van Waer- 
beke et al. 2000; Bacon Refregier & EUis 2000; Wittman 
et al. 2000; Kaiser, Wilson & Luppino 2000). As weak 
lensing maps of the convergence field, k, cover larger areas 
with better signal-to-noise, statistical properties beyond 
just the second moment can be measured. It is therefore 
important to obtain theoretical predictions for the behav- 
ior of the genus statistic for the weak lensing field, which 
we present in this paper. 

The linear-scale fiuctuations of the mass density field 



preserves the Gaussianity of the initial state. However, 
measuring the linear-scale fiuctuations is difficult due to 
sample variance, unless the area covered is extremely 
large. In forthcoming weak lensing surveys, the conver- 
gence field is likely to contain non-linear effects of grav- 
itational clustering. The genus curve is distorted by the 
non-Gaussianity induced by gravitational evolution, which 
precludes a simple comparison of the observed genus with 
analytical predictions. Since this effect is dependent on the 
cosmological model, it is of interest to establish whether 
forthcoming surveys will allow for an adequate measure- 
ment of the non-Gaussian genus curve. 

We show in this paper that the effects of gravitational 
clustering on the genus curve can be measured from high 
quality data. For this purpose, we analyze the ray-tracing 
numerical simulation of weak lensing developed by Jain, 
Seljak & White (2000). The noise due to the finite num- 
ber of galaxies used to probe the weak lensing field reduces 
the signal-to-noise on very small scales. Therefore we have 
to balance sample variance on large scales with discrete- 
ness noise on small scales to find the best length-scales for 
measuring the genus. We study this issue by simulating 
the random noise expected in observational surveys. 

This paper is organized as follows. §2 introduces the 
genus statistic. In §3 the main results of our analysis from 
ray tracing simulations are presented. The genus curves of 
simulated noise-free convergence fields are shown in §3.1 
and those of noisy fields in §3.2. We conclude in §4. 

2. GENUS STATISTICS AND WEAK LENSING 

We consider the local convergence field k and calculate 
genus statistics of this field. In the weak lensing regime, 
K <C 1, an area in the lens plane is magnified by a factor 
I + 2k compared to that in the source plane. The local 
convergence field can be reconstructed from the shear field 
(Kaiser & Squires 1993). For a distant observer, under the 
Born approximation, the convergence along a line of sight 



can be expressed as a projection of the density perturba- 
tion S, 
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where fig is the present value of the mass density parame- 
ter, Hq is the Hubble's constant, o(x) is the scale factor at 
the conformal lookback time x, c^aIx) is the comoving an- 
gular diameter distance, x the comoving radial coordinate. 
If the radial distance is measured in units of Hq^ , we can 
drop the Hq in the above equation. Integration is per- 
formed along the unperturbed light path. In this case, the 
statistics of the convergence field can be computed from 
the statistics of the density field. 

To quantify the topology, we find contour lines of a 
threshold value labelled by i^. For a Gaussian random 
field, the threshold v is identified so that the field takes 
the value k = va, where a is the rms of the field, a^ = ('*^)- 
The genus is defined as: (number of contours surrounding 
regions higher than the threshold value) - (number of con- 
tours surrounding regions lower than the threshold value) 
(Adler 1981; Coles 1988; Melott et al. 1989; Gott et al. 
1990). The genus per unit area G2{i^), as a function of the 
threshold is known as the genus curve. 

The genus curve is analytically given by (Melott et al. 
1989): 
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where (fc^) is the square of the wave number k averaged 
over the smoothed 2D power spectrum. Both sides of the 
above equation have dimensions of [steradian]^^. If the 3D 
power spectrum is given by a power law with spectral in- 
dex n, and a Gaussian window function T4^(r) oc e"'" ' '^^^ ■* 
is adopted, then (fc^) = {n + 2)/(26lf^) for n > -2. 

The above analytic results are not expected to hold for 
ongoing weak lensing surveys because the survey areas are 
not large enough to probe the linear scales. However, it is 
known that the shape of the genus curve of mass density 
fields is not strongly affected by nonlinear gravitational 
evolution provided we use a rescaled labeling of ly (Gott, 
Weinberg & Melott 1987; Weinberg, Gott & Melott 1987; 
Melott, Weinberg & Gott 1988). For a non-Gaussian field, 
the threshold is defined through the fraction of area / on 
the high-density side of the contour lines, which is related 
to the corresponding values of ly for a Gaussian field as: 

f = ^ I e-'^'^dt. (3) 



For a Gaussian field, this definition gives the same result 
as specifying the threshold directly in units of a. For a 
non-Gaussian field, it compensates for the horizontal shift 
of the genus curve due to the shift of the probability dis- 
tribution function. We adopt this prescription for setting 
the threshold values of k in the results shown below. 

3. NUMERICAL RESULTS 

3.1. Genus without noise 

We use maps of convergence fields from the ray-tracing 
simulations of Jain, Seljak & White (2000). The source 
galaxies are taken to be at z = 1. We use two cold dark 
matter models: a "tau" model with Op = 1, and shape 



parameter F — 0.21, and an open model with f^o = 0-3, 
F ~ 0.21. The fields we use are 166.28 arcminutes on a 
side for the tau model, and 235.68 arcminutes for the open 
model, sampled with 2048x2048 grids (see Jain, Seljak, & 
White 2000, for further details). Figure 1 shows a k map 
for the open model. 

The genus curves are calculated using the CON- 
T0UR2D code and related routines, which were kindly 
provided by David Weinberg (see Gott, Melott & Dick- 
inson (1986) and Melott (1990) for details of the meth- 
ods). The key relation used is a 2-dimesional analog of 
the Gauss-Bonnet theorem, which expresses the genus in 
terms of the integral of the curvature along a contour line, 
G2 = J C ds/2TT, where C is given by the inverse radius of 
curvature r^^, and is negative or positive, depending on 
whether a low or high density region is enclosed. For each 
model, genus curves are obtained for 7-10 realizations. The 
error bars on the genus curves shown in Figures 2 and 3 
give the la deviations among those realizations, and repre- 
sent the Icr error expected for surveys of our sample area. 

In the left panels of each plot in Figure 2, the genus 
curves per steradian of the simulated convergence fields 
are plotted for the open and the tau models. Each plot 
has a different smoothing angle, 9{ = 1', 2',3', 10'. In the 
plots, thin curves are the Gaussian predictions with nor- 
malizations fitted by minimizing chi-squares of the simu- 
lation data. The smoothing angle of 10' corresponds to a 
scale on which linear theory can be applied (e.g., see Fig. 
3 of Jain & Seljak 1997). Within the sample variance error 
bars, the genus curves of the simulation data are consistent 
with Gaussian predictions on this linear scale. 

On smaller, nonlinear scales of 9{ = l',2',3', the genus 
curves of the convergence field show a "meatball shift": 
a dominance of the number of isolated clusters over iso- 
lated voids. This shift is consistent with topological 
analyses of the galaxy distributions (Gott et al. 1989, 
1992), and is a common feature of cold dark matter 
models (Springel et al. 1998). To quantify the degree 
of the meatball shift, we define a meatball-shift statis- 
tic, S = J G2{i')diy/ J \G2{i')\di', where the integration 
is limited from v ~ —2 to +2. Likewise, the departure 
from Gaussianity can be quantified by a statistic like: 
5*2 = / \G2{i^) - Ggauss(j^)|// |Ggauss|- We discuss some 
results in the following sub-section. 

The absolute values of the genus number of the high- 
convergence side (v > 0) is larger for the open model, and 
smaller for the tau model. Since our open and tau mod- 
els have the same initial power spectrum, the difference 
in the genus number is due to the difference in the non- 
linear gravitational evolution, and, more importantly, in 
the space-time geometry. A volume factor V{z) for a unit 
solid angle is larger in the open model than in the tau 
model. Nonlinear structures also form earlier in the open 
model. Both these effects cause the topology statistics for 
the open model to have different non-Gaussian features 
from the tau model. 

3.2. Genus with realistic noise 

The weak lensing fields in observational surveys con- 
tain noise due to the finite number of source galaxies with 
randomly oriented intrinsic ellipticities. We build simu- 
lated noisy maps of the convergence field by drawing each 



component of the intrinsic ellipticity of a galaxy from a 
Gaussian with variance (e.g., van Waerbeke 2000) 

2 
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where ct^ is the rms amphtude of the intrinsic eUipticity 
distribution and Ng is the mean number density of source 
galaxies in each cell. In the following, the number density 
of source galaxies is 30 per square arcminute, and their 
rms intrinsic ellipticity is a^ — 0.2. 

Figure 3 shows the genus measured from a k field re- 
constructed from pure Gaussian noise which corresponds 
to a map in which galaxies are randomly distributed with 
random ellipticities without weak Icnsing signals. In fact, 
any white noise has the same genus curve of equation (2) 
regardless of its variance. The agreement with the ana- 
lytic curve for a random Gaussian field (Eq. 2) tests our 
reconstruction scheme and genus measurement. In Fig- 
ure 2 the right panels in each plot show the genus curves 
of simulated signal-|-noise maps of the convergence fields. 

For the large smoothing angle 9{ = 10', the genus curves 
are again consistent with the Gaussian predictions. For 
9{ = 3', the noise effect is small, but the meatball shift 
is reduced from the noise- free maps. For 6f — 2', the 
meatball shift is completely erased by noise except in the 
vicinity of i^ = 0. For Of = 1', there are interesting 
features due to the noise. The relative effect of noise is 
stronger for the open model because the lensing signal is 
small [see Eq. 1] . The high-convergence region of j^ ^ 2 is 
not much affected by the noise. This is because most of 
the high-convergence peaks come from strongly clustered 
regions. Around t^ ~ 1, the noise contribution equals the 
signal. For the low-convergence region of —i/ ^ 1-2 the 
genus is dominated by noise because relatively large ar- 
eas of low-convergence (which corresponds to voids) have 
noise troughs superimposed on them. Such complex situ- 
ations make the genus curve have larger amplitude than 
the Gaussian at negative i' for 6f = I'-the opposite of the 
case for the noise-free maps. There is an interesting coun- 
terpart to this effect in the statistics of peaks studied by 
Jain, van Waerbeke (2000); they also found an increase 
in the number of negative peaks (relative to the Gaussian 
case) in the signal-|-noise maps. 

These effects on the genus curve discussed above can 



be quantified using the statistics giving the meatball shift 
and departure from Gaussianity, S and 82- For the meat- 
ball shift statistic we find that the open model shows a 
significant negative (sponge- like) shift ioi 6 — I': S = 
—0.08 ±0.01. It shows a positive (meatball) shift for large 
angles, e.g. for 9 = i', S = 0.06 ± 0.02. The tau model 
shows a similar trend but with lower amplitude. This is in 
contrast to the noise-free field, for which there is a positive 
meatball shift which decreases in amplitude with increas- 
ing angular scale. 

4. DISCUSSION 

We have analyzed the genus curve of weak Icnsing con- 
vergence fields. On large angular scales, the genus curves 
are consistent with the random Gaussian prediction. On 
scales of a few arcminutes, nonlinear gravitational clus- 
tering produces a non-Gaussian "meatball" shift in the 
genus curve. Its non-Gaussianity also differentiates the 
open and fiat cosmological model for a range of smoothing 
angles and threshold values. The noise dilutes the meat- 
ball shift but it does not change the genus curve to the 
random Gaussian one, even on the smallest scale 9f = 1'. 
Over the range of scales tested, 9{ = I'-IO', the open and 
flat cold dark matter models can be distinguished even in 
the presence of noise. 

We have analysed simulations of survey areas of about 
10 square degrees. When survey areas approach 100 square 
degrees, the error bars shown in Fig. 2 should reduce by 
about a factor of three. Thus, provided systematic errors 
can be kept to below the level of the random noise modeled 
here, the genus statistic can provide a valuable means of 
measuring the non-Gaussian features of the lensing signal 
and distinguishing cosmological models. Further work to 
measure other topological statistics and simulate a wider 
set of cosmological models would be useful. 

We are very grateful to David Weinberg for providing 
us useful routines to calculate the 2D genus curve and for 
stimulating discussions. We thank an anonymous referee 
for suggestions that helped improve the paper. T. M. ac- 
knowledges support from JSPS Postdoctoral Fellowships 
for Research Abroad. B. J. acknowledges support from 
NASA through grants NAG5-9186 and NAG5-9220. 
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Fig. 1. — Comparison of the reconstructed and input maps of k. The left panel shows the input k maps for a field 3 degrees on a side, 
and the right panel shows the same field reconstructed from simulated ellipticity data that includes the noise due to the intrinsic ellipticities 
of galaxies. 
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Fig. 2. — Genus curves per steradian for simulated convergence fields. The four plots use smoothing angles 0f = 1', 2', 3', 10'. In each 
plot, fields without noise (left panels) and with noise due to the intrinsic ellipticitics of galaxies (right panels) are shown. Symbols and thick 
lines correspond to the genus measured from the simulation data. Thin curves correspond to the prediction of random Gaussian fields with 
normalization fitted to the simulation data. Two cosmological models are shown: the open model (solid), and the tau model (dashed). 
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Fig. 3. — The genus curve of the pure Gaussian noise. This curve corresponds to a map in which galaxies are randomly distributed with 
random ellipticities and there is no weak lensing signals. The error bars are given by the Icr deviations from 10 realizations. The area is the 
same as that simulated for the open model and the smoothing angle is 9f = 1'. The solid curve shows the prediction for a random Gaussian 
field. 



